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ABSTRACT 


A  mathematical  model  is  developed  for  the  Class  V  Flextensional 
Underwater  Acoustic  Transducer  Shell. 
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INTRODUCTION 


'’’he  flextensional  underwater  acoustic  transducer  concept  is  pres¬ 
ently  undergoing  extensive  analytical  evaluations.  In  general*  the 
various  flextensional  designs  can  be  placed  in  one  of  five  different 
classes  [ l] .  This  report  describes  an  analytical  model  developed  for 
the  Class  V  flextensional  underwater  acoustic  transducer shell.  A  picture  of 
this  concept  is  shown  in  Figure  1  and  a  detailed  sketch  of  a  typical 
design  is  given  in  Figure  2. 

Flextensional  transducer  designs  of  the  general  type  shown  in  Figure 
1  was  originally  proposed  as  having  possible  applications  as  asonoV-cy 
transducer.  As  shown  in  Figure  2  this  type  of  design  consists  of  two 
shallow  spherical  shells  bonded  at  a  boundary  and  a  thin  piezoelectric 
disk  Joined  at  this  boundary  by  utilizing  an  epoxy  cement.  The  piezo¬ 
electric  disk  is  isolated  electrically  from  the  two  shells  by  removing 
the  silver  electrodes  beyond  the  region  of  contact  between  the  shells 
and  the  ceramic  disk.  Sufficient  epoxy  is  applied  so  as  to  fimly  at¬ 
tach  the  disk  to  the  inside  shell  boundary.  Two  small  holes  are  drilled 
through  the  shells  and  serve  as  entrance  for  the  electrical  leads  to  the 
electrodes  plateu  onto  the  ceramic  disk. 

Although  this  class  of  flextensional  designs  is  designed  primarily 
to  be  used  as  a  sonobuoy  transducer  there  are  special  environments  where 
this  type  of  transducer  could  be  used  mainly  as  a  source  of  acoustic 
energy.  In  this  instance  an  additional  clsmping  load  would  have  to  be 
applied  around  the  boundary. 

If  the  Class  V  type  of  flextensional  design  is  to  function  satis¬ 
factory  as  a  sensor  then  the  flat  portion  of  the  response  curve  needs 
to  be  as  broad  as  possible.  Associated  with  most  attempts  to  increase 


FIG.  1.  PICTURE  OF  FLEXTENSIONAL  SHALLOW  SHELL  SONOBUOY 


the  sensitivity  of  the  shallow  shell  concept  is  a  reduction  in  the  sys¬ 
tems  fundamental  resonant  frequency.  Of  course,  a  reduction  in  the 
fundamental  resonant  frequency  reduces  the  usable  frequency  range  of 
the  concept.  Ac  empirical  equation  has  been  derived  that  effectively 
predicts  the  sensitivity  of  this  concept  below  the  fundamental  resonant 
frequency . 

If  the  total  performance  capacity  for  this  type  of  transducer  design 
is  to  be  fully  realized  then  it  is  necessary  that  a  detailed  mathematical 
model  be  developed.  The  purpose  for  this  report  is  to  present  an  anal¬ 
ytical  model  that  can  predict  the  dynamic  characteristics  for  this  type 
of  sonobuoy  shell  design.  Of  course  once  a  dynamic  model  of  the  shell 
exists,  then  conbining  such  a  model  with  the  solution  developed  for  a 
thin  piezoelectric  dish  with  nn  arbitrary  impendence  on  the  boundary  [3] 
will  result  in  a  math  model  for  the  complete  system  in  air.  If  in  ad¬ 
dition  the  external  acoustic  loads  are  determined  by  utilizing  a  numeri¬ 
cal  technique  such  as  has  been  developed  by  Hess  IM,  then  a  complete 
math  model  'll!  exist. 

The  math  model  describe!  by  the  main  bod;-  of  this  report  assumes 
that  the  edges  of  the  shells  are  horizontally  guidsd-pianed .  In  at¬ 
tempting  to  determine  an  empirical  equation  that  would  consistently  predict 
the  receiving  sensitivity  of  the  shallow  shell  concept  considerable  dif¬ 
ficulty  was  encountered.  One  reason  for  this  difficulty  was  the  inability 
to  establish  the  degree  of  c leaping  between  the  surrounding  shells  and 
the  ceramic  disk.  Possible  causes  of  tbis  variation  ate  the  size  of  the 
shell-ceramic  contact  area,  the  variation  in  the  stiffness  of  the  bond 
joint  which  holds  the  ceramic  (  ■„  A  the  two  shells  together  and  the 
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variation  in  the  thermal  expansions  between  the  shells  and  piezoelectric 
disk  during  the  curing  stage.  It  has  been  initially  assumed  that  the 
bond  joint  acts  more  as  a  pinned  boundary  than  as  a  clamped  one.  Also 
from  a  practical  standpoint  it  has  been  necessary  to  taper  the  edges 
of  the  shell  as  is  shown  in  Figure  2. 
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STRAIN  ENERGY  IN  A  SPHERE 


Static  Potential  Energy  of  a  Sphere 

Using  the  results  obtained  by  Langhaar  [5]  ,  one  can  specialize  them 
to  obtain  the  strain  energy  of  a  sphere  or  spherical  section.  Noting 
Fig.  3  which  is  the  same  as  McDonald's  [  6  ]  ,  one  can  write  the  parametric 
eqs.  of  a  sphere.  These  are 


s  Ck^. 


© 


j 


-  Ck..  0  c  j&'C.  Q  ^ 

and  =•  (X. 

In  eqs.  ;i),  the  parameters  are  choosen 
PC-  0  and 


(I) 


(2) 

Using  Langhaar1 s 


and  the  radius  of  the  sphere  has  been  denoted  as  a 
approach  along  with  (2)  above,  one  finds 

A  ,  U  ^  A  V  M  ^_l.  .at 


f.X-*  V'«>.**X,  +  V* 


and  from  (1) 


so  that 


1  -  (X-  0  og-g  S  . 

0 

=■  (X  i 

fi  . 

0  * 

E  a  Os* 


&  ) 
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Likewise,  one  finds 


where 
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Figure  3 

Coordinate  Systea 


1G 


so  that 


P  '  S  A 

Vtr  -  cv.  C? 


Nov,  neglecting  the  quadratic  terms,  and  the  terms  (e  -  2a‘;/E)x 
-  2 

v  and  (8  -  2g  /G)  v  one  obtains 
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By  looking  at  (6)  and  (7)  it  is  seen  that  e  ^  g  are  needed  so  that 
specializing  his  ellipsoid  of  revolution's  results  (a®b),  one  obtains 

-  CX_  £  _/xjL/yv_  ©  Si  j  ^  -  Ow  H  S' 4^ 

^  s-  CX-  E  ^  c>?y^  0  J 


hence 


_  -1/5. 

Q-  =  -  CXw 

E 

-  l/o3. 

£ 

.ae. 


On  using  13)  and  (5)  equations  (8)  and  ' 9 )  become 

&  .  a  ^  y/\  .  £ 


©  (_y*V~  (Si.  ;  a  -  ©-O-k^vv.  ^ 

;r  -  c-so  0  . 


e.-  -  -<*- 


^  -  -cx_  .AjL/rv_  0  . 
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(12a) 


Hence  using  (3),  (5  ,  and  (13)  in  (6)  gives 
CX,l  f  *1 KT  )j 

^,S  ©*  -Hat  Aa^v  e)^-u  ©, 

VQ^-Cu.^  t  &  -  V  oa-3@)  , 

Likewise  (7)  gives 

***  “©  ‘  (12b) 

* "ur^  a*  w-,.  0-o<voft-  i ^ q  ^ 

=  'Urfi  og*.  ©  , 

£qs.  (12a  J  and  (lab)  tgree  with  (21)  and  (2?)  obtained  by  Langhaar  [5] 
on  pg.  187. 

Therefore  the  energy  due  to  stretching 
for  a  sphere  is  as  follows: 

"ur * i±_  C**  (®f(u.  ■f-vj-')  •+■ 


r  s*( 


+  ©)*  t  aO  tAC  ©Cu_  .  +  7xn(jr^tlA.c»a  © 

+  -ur  ©)+  I  (.i-iJ)c-up“0(.t^-+v-yJ^.B 

--V  ©rj  V\  ©  £e  .  (U) 


It  should  be  noted  that  (13)  agrees  with  eq.  (2.1.1)  of  McDonald's  (6  ) 
paper.  (McDonald  uses  u,  v,  and  w  meaning  u,  v,  w). 

The  bending  energy  ,  V?,  can  be  similarly  obtained,  1+ 
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Eqs.  (  lUj  is  also  seen  to  agree  with  McDonald's  [6  ]  eqs#  (2.1.2). 
The  potential  energy  of  the  external  forces,  ft,  denoted  in  Fig.  1  by 
y.y,  i  readily  give 

«T  © 


J  L*  -of 


J  s' u  ,A_ 


+  4  vt  2  “W^  -£x^v.  &  (JL0r  c^_(^ 

(15) 


Hence,  using  (13),  (lM  and  (15),  the  total  potential  energy^ 


of  the  system  is 


\/»  Or  +  Or  -f  -A- 

I  A 


(16) 


Exact  Solution  In  Circumferential  Direction 


The  following  middle  surface  displacement  fields  are  assumed: 

CO 


uJ.9  £)  =  u»  U)  r  'vr\.  Q 

*  1**0,., a,...  -*«  •> 

irC&jQii  2L  ir  (. 6) xs. Si  . 


and 

-u r(8i  5.  1>U"  C  0^0  ITS.  &  . 

.  ^nrs- 

1U 


(IT) 


These  expressions  can  then  he  substituted  into  equations  (13)  and  (±4j. 

To  use  (lT)in  (13  ,  and  (lM,  one  should  note  the  derviatives  with  respect 
to  0,<(>  of  (IT).  To  differentiate  the  series  in  (IT)  one  assumes  all  the 
necessary  conditions  as  noted  in  Widder's  [  T]  Advanced  Calculus,  pg.  305, 
extended  to  two  variables. 

Hence,  c0 

rS~& t  O  -  IA_  L0)  CJS'Ci  irv\.lQ 
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nr  L©  Q)  -  SL  -vyn.  V  Lb)  r  TTV  Q 
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•"Ws*  '  > 


1>J  L  &)  s.  _Al^.L0)  ,  . 
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(18) 


and  2.-JYVU7  L©)^ _ .  / 

*r/N.  >  / 


Consider  the  first  term  of  (13)  and  using  (  T ) and  (l8), 

(u-  4-  'Iat)^'  a  r  Sl  C la  '  +^xr 

d  L  ^  ^  J 
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Taking  the  integral  w.  r.  t.  $  inside  the  ,  sign  and  utilizing 
orthogonality  gives 

C  Upw)  «ia.  S _  u' 


0,1,8,... 


Cu.'  +V4)t. 


1A.  ivv. 


Similarily  one  can  work  with  the  remaining  term  of  equation  (13). 


Therefore  the  strain  energy  due  to  membrane  becomes 

.  .  „  ,6  co  r  ,  . 


l-V  J 


^ C  5:  fCuJ  -v'ur 

|*i;  ^  *^*0,1,31,  L  ^  ^ 

■+  t/*  0  -v  -ur  ^^0fV  anJc^©luw‘ 

^nv.  >yvv-< 

+  nv^J  Uv.  v^+  B  +  'UT^-Al^  0  V 
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**  -Yy\.  -WV.  -VA.  _J 

K  Al.wv  0  cSl.0  «  (20) 

Similarly  the  strain  energy  due  to  bending  is 

\S  .5  L ( "*£?+  c  e  s^8 

*0  ^naO,l;A,..- 

Kiat' 

-rrv»  *TV\» 

-  /WvS'  'W1ft)+  s(l-l))  oaC*  8(J*>L'W^Let  fir-  -TA_"W  1  ^ 

>Cl  A../^  &<$-©  .  (21  ) 


A  |S((|*1) 


Consider  the  following  substitutions  for  the  load  terms: 


t 


where  8m  =  Fourier  coefficient  depending  on  the  circumferential  dis¬ 
tribution  of  the  load.  Using  (i7)and  (23)gives  using  orthogonality 


J\.  --TT 


CO 


CX- 


* 


X0  ^=0,),* 


U_  +  f 

>YA»  .  1 

(23) 


Physical  Interpretation  of  Assumed  Displacements  and  Loads 
Consider  the  case  n=0  of  eqs.  (  )  and  (23j.  For  ;17) 

u_ (_©,(£)  =.  u- IS) 
or  (.9 -  0  j 

“d  -W  (.  0  ,  Ql  a  C.&)  •  <2l>) 

From  (2U;  the  displacements  are  independent  of  d>  and  are  thus  rotionally 
symmetrical  about  the  axis  of  revolution.  In  this  case  the  displacement 
of  any  point  is  the  same  as  every  other  point  on  the  same  latitude. 

Thus  one  can  show  a  w-displacement  on  a  cross-section  of  the  spherical 
cap  to  be  something  similar  to  Fig.  •  which  could  be  rotated  about  z  to 
give  the  shell  configuration. 

Likewise,  the  load,  u’3)  reduces  to 

X  *  *)j 

H  -  O  (25) 

and  2  S  ^  *C&)  , 

where  J^'Q)  to fc(6)  are  forces  per  unit  area.  Hence  from  {IV  in  the 
case  of  rotational  symmetry,  there  is  zero  load  in  the  y-direction.  And 
in  a  cross-section  view  which  again  could  be  rotated  about  z,  one  notes 
that  the  loading  is  similar  to  that  shown  in  Fig. 5  . 

One  can  next  observe  the  m»l  term  and  its  effect  on  displacements. 

Eq.  (if)  gives 

17 


U.C  U_  jC©') 

V  C  (Q)  ^  V’  ,  (0)  -Aa^VN^Q;  (26) 

*tlT  C  ©;  &  )  -  tu*;  C©)  C-©^  Si  • 

Thus  if  one  assumes  some  6^,  where  O<0,<§,  to  be  the  latitude  to  be 
observed,  and  observing  only  the  v-displacement ,  then  the  v-displacement 
as  a  function  of  $  on  the  0*0^  latitude  is  as  shown  in  Fig. 6. 

It  should  also  be  noted  that  eqs.  (26)  relates  the  stretching  type  of 
motion  in  some  manner.  Thus  when  u  and  w  are  taking  on  their  maximum 
and  minimum  values,  the  v-displacements  (circumferential)  is  zero; 
however,  v  is  maximum  when  u  and  w  are  zero  so  as  to  allow  u  and  w 
to  take  on  their  maximum  and  minimum  values  by  stretching  at  this  point. 

One  can  likewise  carry  on  this  analysis  and  look  at  other  values  of 
m,  however,  it  is  felt  that  a  good  physical  grasp  should  have  already 
been  obtained. 
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Figure  6 


.ypical  v-displacemem.  for  and  some  6«-d 


APPROXIMATIONS  FOR  THE  DERIVATIVES  IS  THE  BOUNDARY  SEGMENTS 


The  crown  point  of  the  dome,  the  point  of  zero  meridian  (co-latitude), 
is  considered  a  boundary.  Hence  a  statement  must  be  made  about  the  dis¬ 
placement  at  this  node. 

Two  cases  are  distinguished.  The  first  is  that  in  which  the  vi¬ 
brations  are  rotationally  symmetrical  about  the  axis  of  revolution.  In 
this  case  the  displacement  of  any  point  is  the  same  as  ever;'  other  point 
on  the  same  meridian.  This  case  arises  when  m  ■  0.  Therefore,  the  only 
possible  motion  at  the  crown  is  one  in  which  the  crown  node  is  displaced 
in  a  radial  direction  only,  with  no  accompanying  tangential  movement,  and 
the  slope  remains  zero.  This  is  expressed  mathematically  as 


m(DV  ir(0'U  O 

\  i 


(27) 


McDonald  [  6}  chooses  to  use  the  finite  difference  expression  for 

2 

the  first  derivative  having  an  order  of  error  of  h  while  the  second 
derivative  expression  is  the  one  with  an  order  of  error  of  h.  Since 
he  starts  at  *  •  0  and  proceeds  positive  along  the  segments  then  the 
forward  difference  is  thus  the  one  used.  Hence  for  the  general  dis¬ 
placement,  the  first  and  second  derivatives  are 


2: 


(28) 
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o  H 


and 
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^WW' 
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(29) 


vhere  h  ■  208  »  2/a.  Hence  using  the  B.  C.  given  by  (Zftin  connection 
with  (26 .'and  (29), 


u_  =  1/  a  o  . 
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(30) 


so  that 
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Also  since  w  ' 
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(32) 


so  that 
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(33) 


In  the  second  case,  since  the  displacements  are  in  general  a  function  of 
$,  the  longitudinal  angle,  then  they  must  be  zero  at  the  crown  or  the 
displacement  at  the  crow  would  be  different  for  each  value  of  Hence 


-  -v  C  -  w*  Lo)  -  0  j 

so  that  eqs.  (2)  and  (3)  become 


The  boundary  conditions  at  the  lower  edge,  9,  depend  on  the  type  of 
support.  In  the  case  of  the  clamped-edge  dome  studied  by  McDonald  [6  ], 
the  3  displacement  components  and  the  derivative  of  the  radial  displace¬ 
ment  must  vanish,  i.e. 


u.  (s')  i  -v  C  D)=  -w(  &)  -D,  \ 


and 


d  B 


(36) 


e=e 


Using  the  seme  order  of  errors  as  in  the  crown  boundary  condition, 


the  backward  finite  difference  expressions  are 
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H  A)-a  AJ-|  * 

vll(0)=i^(1/  -a-ir  ). 

H  ^-2-  /\j- 1 


luJ  (  0 )  -__L  “UT 

H  AH 


%xr  )»0 
A/- a  ' 


so  that 


■ur"(6)* 


i  * 

-1.  <=*  -U-r 

A  M' 


One  possible  boundary  condition  for  our  problem  is  that  of  a  radial 
guided  boundary  with  the  pinned  condition( there  might  be  some  torsional 
spring  like  effect  but  for  now  it  is  neglected).  This  arises  from  the 
physical  fact  that  two  of  the  spherical  caps  are  placed  back-to-back  so 
that  it  is  assumed  that  motion  in  one  is  the  same  as  motion  in  the  other. 
The  B.  C.  is  shown  more  clearly  in  Figure  7  which  denotes  a  cross-section. 
Thus  from  the  physical  problem,  the  case  of  rotational  symmetry  is  the  only 
one  of  interest,  hence  it  is  assumed  that  m  =  0.  In  this  case  the  boundary 
conditions  on  the  lower  edge  are 


moment  at 


rnj 

© 


(1.0) 


and 

displacement  in  z-direction 

6m  § 


(1.1) 


There  is  no  boundary  condition  on  v  because  v(6,  $)  =  0. 

The  B.  C.  given  by  (h 0  /  will  be  discussed  later  on;  when  ex¬ 
pressions  for  moments  are  written;  however,  the  B.  C.  given  by  (1*1)  can  be 
considered  here. 

With  respect  to  Figure  8 


displacement  in  z-direction 


-  -w  c&a  0-  t/L  -W  & 

Ai  Aj  > 


0^6 


(1*2) 


and 

displacement  parallel  to  xy  plane 
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Displacements  at  guided-pinned  edge 


2  6 


Using  (1*1)  with  (1*2)  results  in  one  mathematical  boundary  condition. 


-ur  »  u.  "ttun  © 
A/  /V 


(MO 


or 


(1*5) 


Another  possible  boundary  condition  is  that  of  the  radially  guided 
clamped  end.  This  condition  requires  that 


displacement  in  z-direction 


e--5 


and 


5>-uj- 

$>  0 


(1*6) 


The  first  condition  again  gives  rise  to  either  (MO  or  (1*5).  The  second 
condition  by  use  of  (37)  gives 


W  >  H  -1A7 
N  3^  N-l 


~  1  -m-T 

“T  N-a^ 


(1*7) 


Physically,  the  boundary  condition  is  somewhere  between  the  clamped 
and  pinned  cases  of  the  guided  endsj  however,  how  much  between  is  unknown 
since  the  exact  effects  of  the  epoxy  in  terms  of  an  equivalent  torsional 
spring  effect  is  unknown. 
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STRESS  RESULTANTS 


The  stress  resultants  and  hence  the  stresses  in  the  shell  can  be 
determined  from  the  displacements.  This  is  done  by  the  fact  that  dis¬ 
placements  in  the  shell  are  known  and  have  been  shown  to  be  linear  in  z. 
From  displacements  one  can  determine  strains,  and  then  by  use  of  Hook's 
law  for  an  isotropic  media,  one  can  determine  stresses  as  a  function  of 
displacement.  The  stress  resultants  are  then  found  by  integrating  the 
stresses  across  the  thickness.  McDonald  [6  ],  who  references  Vlasov  [8  ], 
gives  the  following  eqs.  for  the  stress  resultants: 

q  =  D  wBe+' L  a  ^  B+1^  *3  > 

ks)+  Cu-Cj&t  &+  iTq  ca*  &+vr)j } 


(1*8) 


where  M_.=moment  per  unit  length  acting  in  the  plane  of  a  latitude,  i.e. 
in  the  direction  6-direction, 

M  =moment  per  unit  length  acting  in  ^-direction, 

w 

Kp*aid-plane  force  per  unit  length  acting  in  0-direction, 

N  »mid-plane  force  per  unit  length  acting  in  ^-direction. 


Eh 


3 


12(i-\T) 


and  K= 


Eh 


i(i-vc: 


The  positive  direction  of  the  stress  resultants  is  shown  in  Fig.  9. 

How  substitution  of  displacements  in  their  assumed  series  form, 
eqs.  (  17)  into  i*8),  gives 
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N\  =•  o  %.  [V 
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ce)  cA.  a) 
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N\  '-ol  f  -0-ur "(©)+('  ^<s^J.&)cuc.&8+-  J 
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+  -u/He)  of*x.  b  j]  goo  j  r9) 
N [\j. (eVWeVOi^.  Ce)e>rft0  \ 

6  -ta»0  L  ■*'-  o 

+  v^e)  c^^Le)tv^(e)jJc^'w^J 
n  »*  iL  f-Olu.1  (j&)+'wL^e))+Cu-  (e)c®t  e- 

/vas.0  L  >*yv.  *^r>“  ^  n 


/vAs.0  l_  >*»v.  'TfV- 

C©))  ce< 

—*s.  -WV. 

’•or  a  spherical  cap  experiencing  symmetry  deformation,  one  considers  on 
the  m=C  term  so  that  from  (t^) 


,09°V 


N\ e-  or->r "  (9)tt)v Xe)  oefc.  0^]  , 

M\  =0^uj-"(e))+-W,(©)cjet  ©]  ,  < 

N  -  is  ^  uJ  (.©Yt  -w  Is) + t)  L  ^.(.e'lceC  ©+  w  (e))J ; 
dl\|  ^  s<.  jjoCuJCe'li  ur(.9)+  i^ie)  c^C8+*  w-  Wl ,  / 

ere  P>  -  g^_  .  K  _  jEX  1 

”  Ito^v*  * 

Ihe  stress  resultant  of  the  iwme  point  can  be  vritter.  from  (“S)vhiie 

.  (»"*■,  the  l our. i ary  conditions  at  the  dome.  To  that 

%  «  Dw“(0), 

°  StO 

Mk.  I  zO0ur"(e)  . 

11  Q*-0  ' 


noting  e^s 


(52) 


N 


0 


0  s.  o 


-H  L©)  i"  ^(©yj  ; 


and. 


|s  -K  [0  u)Ce')+-Liri)y"(ve^ 


Using  the  finite  differences  approximations  and  boundary  conditions, 
eqs.  (51)  and  (52)  become 


^V0 

LN\)  .  -a-ur  f  'Ur  ) 

6  g-0  *  0  '  •* 


(53) 


and 


^%Va° 


D. 

D- 


(.Ta/.  -  ^W,+  T/fi)  * 


(5U) 


For  an  interior  ccint,  the  stress  resultants  or.  the  kth  segment  are 


given  by  (b8)md  to  be 

<r 


(.^  -q-  p 

+=^U+-x)c^-^-WK+l  } 

(55) 

jjO-^  Wj^- 

<r.  ,1, 


(NO. 

<X  ^  -  k  ' 


and 

How 

different 
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This  pro 


A#f> 


(NJ  • 


.  K 
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&  tfK 


ina  k--i 

4i  <j-;K  *• 


U.  +"i«- 

Ic,  -k-n 


+  J  ( 56 ) 

r  x)°s)u,^4-  *1^  u. 

*»  I  v  v  <— » 


yK'H 


4-1  !+*i))ixr^ 
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for  the  lower  edge  four.dary  condition  the  following  finite 
;e  approx  inuit  ions  are  utilized. 


%'n  rU  J 

+  loN'l  ^  * 


redure  yields  finally  the  following  set  of  equations: 

(sNv  ')  .  -  ^r~  f  <~£^-  ©)  -txT 


b  a*N 


-  a  4--  ■> 1  = 

(  n\  )  „  -S-  M>-a  cs^eW 

<s^*N  ^  L  x  N-,  A  j 


\  -  1  a 

l0^c-cj*St-6  )w  +0«>  t*r  , 
N  Ni-l 


CN  Y 


^  L 


=*.  \A. 
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Energy  Expressions  for  a  General  Interior  Point 
Lagrange's  equation  nay  be  written  as 

_aL  (-SL-  ,^)  -  ,  (to) 

where  L  =  T  -  V  is  the  Lagrangian  equal  to  the  difference  between  the 
kinite  and  potential  energy  and  q ^denotes  one  of  the  generalized  co¬ 
ordinates  q  ,  q2, . ,  u.. .  Hence  for  a  stationary  system,  T*0  and  V=V 

(q  ,...,qjj,x,y,z)  so  that  (t'j)  reduces  to 

q»V  =  O  .  (6D 

where  q^  denotes  qg,  q^,..^,qjj  or  (N+l)  quantities  and  q  is  the  general 
displacement  denoting  u,v  and  -w-  or  3(H+l)  iisplacements  in  general.  The 
potential  energy  denoted  in  (61) is  the  total  potential  energy,  and  thus 
is  the  sum  of  the  potential  energy  of  each  segment  of  the  dome.  Thus, 

^  Nv.  •  (62) 

dr -0  ® 

Noting  the  types  of  segments  formed  in  Figure  3  one  can  write  (5)  as 

v=  | rD[Lv'VCv^~A'^-  (63> 

The  potential  "nergy  in  each  segment  is  formed  by  integrating  eqs.(20), 
(2i),  and  U3J  for  the  interval  of  each  segment  and  using  the  appropriate 
difference  approximations  for  the  derivatives  as  outlined. 

Since  we  are  interested  in  rotational  symmetry  only,  then  one  can 
simplify  the  above  equations.  Hence  considering  only  m*0  and  V=0 ,  then 
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V  --—L 
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V\.  po  7T 
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^  ( u~  +  VJ  J  -+■  oo*C,  0  (u^ B ^ J 

+  S'u  bC 

-4^  BA 9  (6U) 


<*  l«fc(l -y)tC 
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~ -  TTc*.: 


X'W’  -^.c^w  &(j&c 

8  [ 


•  '.  •re  r.  's  assures  tu  le  constant .  't>-'' 

Starting  v ;  *  h  the  i„r.e  point  and  working  toward  the  second  boundary, 
the  firvard  and  central  difference  apprcxisations  of  0(h)  vi  i  be  utilized. 


In  (68)there  are  three  types  of  integrals.  They  are  defined  as  the  following: 

I.a  BA  9  t 

\  )  -  ) 


y 


d#- 


'/*=* 

X  a-  ft  A  ^  j 


(69) 


y 
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-  r  ^ 
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O’ 
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Evaluating  the  integrals  of  (69)  yields 


X  —  ^ 


1/c* 


and 


X  ft.  s~i^K  sr\^  1/  tTx 

^  *  ft*  1/ 

X  m-Q\  ft  •  I  /ok  “h  Ovy\^"t(X/A^  — L- _ — 
a  - 1/^.  A 


(TO) 


—  Q-/-W —  Lst/VV 


Using  (TO)  *n  connection  with  (68), one  can  thus  obtain  the  membrane 
energy  in  the  jth  segment  to  be 

(xrj .  --  Kia-7t  j Cu3  -a,  u.  u--  +  u3)+(hi>)(°ma.  v- 

'  #  i--o  (  Lh  a--'1  <y'  d  &  d+l  0 


'  LA-.  "W*  t-  r*  W; 
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1  °\)  0»  /<=<+  ^0+”^  u—  •  ^ 

i  J  d  L  d-  d 


}'4 '"'d'J  ”  d" . L' . y  d 

+  "d  ^  w_.  u_,-b  •»  u~a  cxfva  B  •  >V\_.  1/4*.  +  u_  -  [-3l  /u/a6<  ivriH. 

<4*1  v  d-l  dL  d 

+  Uv-teuw.  -w-w^r  *“■ 


(Tl) 


In  a  similar  manner  the  bending  energy  for  the  jth  segment  was  determined 
and  is  as  follows: 
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,  <9 

-l-SL/rv-U. 


J+I  ^ 

^2  -  JU^-tau^  ^£l)  - 


Observe  that  the  membrane  and  bending  energy  expressions  give 
the  stiffness  matrix  coefficients  while  the  external  load  energy  gives 
the  mass  matrix  by  use  of  D'Alembert's  principle.  Hence  one  can  now  obtain 
the  stiffness  matrix  coefficients  of  all  points,  except  the  9  edge  points 
and  those  next  to  them,  by  use  of  eqs.  (71 ),  and  (72). 

Now  apply  D'Alembert's  principle  to  obtain  the  equations  of  motion. 
Since  by  D'Alembert's  principle 

X  “  -  P  k.  ia-  j 

^  ^  ^  ;  (73) 

£  -  -  y  k  -ur  ^ 


and  assuming  steady  state  conditions. 


gives 


and 


^  =  Pk-ur^v 

i  -  o  k.-ura*'"e-<^t . 

j 


(7U) 


Substitution  of  ilk)  into  (15)  gives  for  rotationally  symmetrical  vibrations. 


A.--  I 


oU  /i4  uJ^  +  w  )  6  ck  Bj  Q.  ^ 
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(75) 


'4 


1  wt. 

Now  substitution  of  u=ue  ,  etc.  into  the  membrane  and  bending  energies 
1  W"t 

drops  the  eJ  from  all  expressions.  Hence  (75)  can  be  written  now  for 
the  Jth  segment  as 
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Energy  Expressions  for  the  Dome  Point 

One  can  observe  that  the  two  equations  of  motion  for  the  dome  point 


can  be  found  by  evaluating 


it-  o  . 

o 

First,  one  considers  eq.  (75 )  where  the  only  energies  containing  u^ 

are  (V  )  "f  (fi)  so  that 
J=°  J=0 

VCu..>(V,V..„+  Ul),  _  „  ,  ;70) 


T  V  '  V=.  O  )  !79) 


v  O 

In  equation  ( 78 )  observe  that 


In  equation  (jQ)  observe  that  /  \  /  \  \ 

VUr>(U  ,^0+  CVa^D+  C  A^-.0+LU^  „ 

Therefore  substitution  of  (81)  into  {78)  yields  the  remaining  eqs.  of 


motion  for  the  dome  point  given  by 


t  . 

The  boundary  conditions  to  be  utilized  at  the  dome  point  are  as  follows: 


u lo)  -  xr  to)  =  0 

c)  I  II I  , 


lo>o  ,  j 

McDonald  [  6]  adds  to  this  set  the  following  boundary  condition, 

3©**  ‘©c  0 -1+~r>'l(.0)  5.  0 

This  is  adopted  in  order  to  obtain  a  finite  solution  at  the  pole. 
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Now  rewrite  the  membrane  energy  of  the  dome  point  segment  for  the  in¬ 


terval  O<0<l/a  so  that 
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tv  Y  b-.tr  TT  \  °|r 

0  I  -  'd  Jo  *i-  ^  __ 

4-  U r  /1>C. 


+i/r  )  -iu/rv.  f)  cL©  . 


(85) 


0  jj 

Using  the  boundary  condition  for  rotationally  symmetric  vibrations,  i.e. 
Uq=0,  one  obtains  the  membrane  energy  for  the  dome  segment  in  this  problem 
to  be 

C 'ir  tt  f  ^  i+ 

■  -  I  '  /r  —  i  J  —  — T-" 1  v  it  I  Ls 


:86) 


^'-u  T=vT{ 

4  Ov(  It  d)^,  l/L(  )l  j  I-  j 

Likewise  for  the  bending  energy, 

(V  .  •  ■=  rrk2  L  -uj~r^  -  -ur^w ,f  -w 

~  laTT^r*  Lib-  ~  ~  ^  ' 


l—Mp  ^ 
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(87) 


In  a  similar  manner  the  inertial  energy  for  the  dome  point  results  to 
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p\^o3 


(88) 

Equation  of  Motion  for  the  Dome  Point 

To  find  the  eqs.  of  motion  of  the  dome  point  w.r.t.  w^,  one  uses 

(82).  To  do  this  the  expression  for  lU  )  is  needed.  To  obtain  (U  })  . , 

£  J-i  -  jBi 

it  is  now  possible  to  make  use  of  the  general  energy  expressions  for  a  Jth 

segment  and  specialize  it  for  the  j=l  point.  Before  doing  this,  it  seems 

best  to  note  now  that  a  meridian  line  is  one  starting 

at  the  dome  point  and  with  increasing  6  arrives  at  5,  while  a  longitude  is 
for  constant  6  and  varying  $.  In  Fig.  (10  )  one  notes  that 
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Figure  10 

Cross-section  view  shoving  segment  divisions  and  their  nodes 


(89) 


Thus  using  6^2/a,  one  obtains 
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Hence  v  the  eq.of  motion  of  the  dome  point  is  as  follows: 
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The  final  stiffness  coefficients  are  defined  similar  to  McDonald's  (9) 

k/  / 

i.e.,  S7T  where  ,  refer  to  the  displacement  component  by  which  the  coefficient 

1J  J 

W 

is  multiplied  and  *  designates  the  particular  eqs.  from  which  the  coefficient 
comes.  Hence  (91  )  can  be  written  as 
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Equation  of  Motion  for  J»1  Point 

The  eqs.  of  motion  for  the  point  next  to  the  dome  point  will  now  be 
obtained  since  the  eqs.  v.r.t.  the  ^  and  displacements  cannot  be 
obtained  from  a  general  expressed.  This  is  due  to  the  fact  that  the 
energy  expressions  for  J«0  contain  u^  and  w^ 

The  eqs.  w.r.t.  u^  will  be  considered  first.  The  energy  which  is  a 
function  of  u^  is 

vU.V  Cv.y.o^cv^+CA)-.,  , 


v  93  5 


so  that  the  eqs.  of  motion  is  given  by 

a^|  C»«_,  3*.,  aiA-'  (9U) 

Substitution  of  U  jso»l  ^1-1  e<lua^^on  (91* )  and  evaluating 


the  partial  derivatives  yields, 
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Equation  ( -f  )  is  the  u^  eq  .  of  motion  of  the  point  J*l. 

To  obtain  the  2nd  eq.  of  motion  for  the  point  J»l,  one  obtains  the 
energy  expressions  which  are  functions  of  w^.  Hence 


v(*r,)4r  1)<s, +C\r ^i0ra/^.^cv ^iat(  iV . 


so  that 


p-v  -  5(V si  t  c^df ali  '^-‘r  c^(  V!!  - 1  ■+- cj)tU  =. 

^  oW,  S^oT,  c)  w.  ^ 

S  .  _  x.  ®  * 


+  cILlLL,*'  =0 


<3  tA/J  ^ 


(97) 


In  a  similar  manner  the  equation  of  motion  of  the  J*1  point  with  respect 
to  was  determined  to  be, 

fer  rk  ftsr«  k-i-csr T)*s+(srrV; 


*h(  i 


' r  CST -UT »£, 


where 


-iaT 

i  I 


•V.  7T  ^ 

ST^STc  /^t)  jo<o/  I  tl))  iU-^w  ”>4*.  &/TV  J 

Sr ?5S  ~ TbJhT)^ j^7 6-  C&'iJ^V^T 

+*&<**&*  A—  J4cj; 

S T7* T^-  /-7j  fo/- jyt  U-<J&>Vr* ) 

'J  *■  VA.  ]/c\+~t7  -Ojl /i«  vv.  X< 

^  °<  cJX.  ^V:/^4x-4a^ 

+■  oc  ~  Q^1it^/^c*<)]  ;  (98) 

C  W  <U>"  . .■  ^  ■  w 

-  /  a  -  i^t-v^s w-c&>&)+-5:  a 

"t  "*}  V_  */^S  ^  C-6e£>  A*wA_ 


-t-  trf'w'CooN^  l^o  -  y^o^J  ) 


S  /  3 


•yv/  ^  ^ 

^ 1  >  . 


!i'  v  to  obtain  the  eq»  of  motion,  one  forms  the  energies  containing 


the  factor  v^.  Therefore, 


-V  ±CU  «*V'a 

*+•  £  •  (99 


Thus  the  v„,  eqs.  of  motion  is  given  by 

Again  substituting  for  the  various  energy  expressions  yielas  the  1  ' 

following  v  eq»  of  motion  for  the  J»2  point. 

t ;  u.A+ cs^  $-)  ^3+^  ?w„+ st  r )  ^ 
-rCS^  r")^ra+  ( T')^%+Csr  ^ M 


a  a 


)-U7n  -  0 


where  t  t- l5  A-o»uX<Jj 

S  «*  ^  |j2o*I,  If1>)./>>W-  \Zo*kTj  j  (loi) 

S^o"  5.  Ax  vv  VA  cM&/Uwv^J 


+^C®c^.  VJUA.  ^+l£|sS? 

CiA  W 

-w_ /^cx" 


^  vS  hjZZ 
o  ^  &  -  1-0 


sr  t*  [gC/t^W  [t-V  «»& ) 

oC^  H 

"t  H  /^.  Ai  ij/w.  *y  >&A/fl.'^ 

-ft/  CK  >6-L/A_.)4><“f  <^7)^x  C^'Cj  Z>W\,  i^s-f  -5>,  >5^w^ 

+  DaV  ti3./A  ^«<x  ”.QyY\££lAV  )£}  «•*)"!■  "^T  ^  t  0>»\tftM  ^o*v 

'_J; 

3  JJ 

-l^ttK  r«^  .  H  ,  .  4  ;  o._  ,  6/  1  t  / 

^  25  3  ■=-  /j?(/ -V)oJ?-  I  ^u. yi\.  /o^  /^xs*^/t?'^T  ^  UU-s>\-  Au\.  /&<• 

.  3 ^ */.,  |>  ot.^  •  |  /  <3^ /  '  H /  .1/ 

T  P  OS^.  -<L>k  ,v'v  -f  ^  (.-£►&  /^<  !WvvJ/<?<  -f-  ^  ^  sfrAA*  £, 

+  -  (U^rbuv  y^c*}) 

^■ur  -W  p  ^  ,  (/>/  .  1/  7j*L 3  V  •  \z  1 

-3  A  H  -/aU-l J)aT  \yr  h>,-r^  /^  <S& /oA  Ajl^/^J  ^ 

*—  _  -O  -V.  ij  .  *  \  J 


and  'hX.y^  y~ C  ^  7TZ>K  U^  "  /WW^  /^c  V^<C  • 

; 


It  should  be  noted  that  the  u„  equation  of  motion  can  he  obtained  from 

£ 

the  general  interior  equations  sine*'  ’t  is  not  affected  by  the  boundary  point. 


Equations  of  Motion  for  the  General  j=k  Interior  :  infs 

How  the  eqs.  of  notion  for  a  general  interior  point  are  needed.  These  vil] 
be  obtained  by  considering  the  partial  derivatives 


(102) 


5^-  -  O  (103) 

.  a-v/K 

For (102  the  energies  containing  the  factor  u  are  the  only  ones  of  interest; 


hence , 


v (u-^)  -  (if  y^K-i  1- (M  i)^*  +  ^ -K 
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Thus,  substituting 


■.  i  0. 


y-e-  '*S. 


-r  r)jy  il  UK  +  C* 

oT  6  0)u_  ^ 

s  ^  <, 


\aK  -  0 


uup ) 


The  energies  containing  the  factor  v.  are  as  loluovs: 

t  X 


U  shr.K-1  +(.1/ 


pCUTV  -k,  f  ~.'JU  "vtUTl^-K  t 

'U/*’  ^ 


(106) 

UK-hl 


i/^  v_^r  t/w  ^ 

S  S  (107) 

■+*  ^C-AY^-^-  0. 

=s - 0 

C>  2a/  ^ 

Tcv  tv  obtain  one  c:  the  general  equations  of  notion  for  the  ti=k  point, 
the  various  energies  define!  in  equation  (iCfc  are  substituted  into  eoaution 


t  ■•:  r  - 


equation  o:  notion  for  the  j=k  point  as  follows : 


CS  )u~  +  Ls£Ol^  +C^kTv^i)  Ul  +  Vu/l 

fs  SM  K-|  »Si-|  IS- 


■e(.S*D  w-  -(rrYO *)  -  0 


W 


vh 


here 


Si>-  u.  ob 

^  k- 1  i  ~'T~TT~  -  <=>«. 


A  ~  ^  J 

^  u,  U.  kv^T-rr  _v  •  ,  \/  .  A  •  aK-  *  |y 

*  K  -  /--_|"  )-*■  ■'W/Y.  j"X^\  **  t  Ai>w.  Os.  >o-< 

“  H  1/=*  ZA4  ^  */o<  "t  Aovv^^  -As^rs  V^s. 

tUw'#-Utu. 


^  Km-TT 
o  K  tS-*-l  -  i -nJ 

-  n- W 


~ 

Ok 


ot 


(108) 


-a*.  •  |y  v  '  wl 


■  iDS  I 

oC  -J 
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o/ k-0  ,  afe-L  .  ak^_\ 

"  J?  i^4  <WtiX/V\-  Si  CK  - WuX/fL  a  C*  J 


-+  Os,  ^Kc>xJsrt  o<"  -4_^n.  )4>  +  k/k."t£X>V\.  ^-^T"  *•  k 


&  .  .L 

<=»<■  ^  ^  ,  <=*  *^ 
i*  S<  c3>^ t  o> 


3  'i  I  \ 

^  -v^r  w  {a*>  it  ^  r  c*  .  a^.lj*\  v 1  j 

->  M,  kt>  a  =•/;*//-  4)^  )  /^  .Ujf^  ^  ,/>c-v\ 


4  ^ 

•  “Sb±!l  .  W  !^si  afltti)  '/] 

CX  •+■  ^  06^  c»<  JXKj^  IgA  J  ) 


(TV  k  ik  i<^  l  o  KcxT  u)  .ilu^  <=*-  J~<s<\  /^ 


c^  /  T  <M  «Jf  ,„'/.  -i%  ^ 

O  ,  2  I  — j  0S  i*' 

I^W>1  li(l-J)*r  L  * 

?  U  v?  ^  |c  4  +  ^  J?  ?  A  5  ~iL* 1  )  SM 

•rf  -4*  £«s.  TJ-  :k'  4  "7  x 


ZM 


+-  4-1,  £■£.  ^  / 


,  t  I  C  i-  t  1  .  L.  fa-v)  . — 
]/  f  Xv»  f:  v\ - -H  *-* 
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Energy  Expression  for  the  Lover  Edge  Half-Segment 


The  energy  expressions  for  the  J=N  half-segment  are 

^  c*  (no) 

o' N  ^ 

Ax'aB')  ■+  c^r  BCLL'n+  WN)(jlfic&&4’t*r/t  9(LQ} 


v,y  -T^j  riw-‘Nt*  +c.^<a&(-ii/sjc^0 


«'»v,-Srau'S®. 

***\  j  ^  v  1 


1  ySi 

J  0  0^6>) 


,  _  >2,  ^ 

t5»L' 


8  C  w«XV  />  /yv^  o<SvC>  £»]  &  A- 

.TTqVncJ  U)  Cx  U.^  'V  VJ  pg  ~\  A^v\-.0£^_fa  * 

f  J .  L  mu 


IJU  u)  C  u.^ ^ j 

^  ^b-)U L  (ii 

Fbr  the  lower  edgejthe  finite  difference  approximations  to  be  used  are 

.  CX  /  v 


%  "  * 


(112) 


^^VVWi>' 

Substituting  the  finite  difference  approximations  and  evaluating  yields, 

£  J  [*“-«  -  SU-N“-N 

+  <*-L\+v)LLLnwu-  U.  (113) 

+  [a(i+o)U.NurN  ^  o-<(.u.*-  ixNu.^]c^S 


+  IL 


a  o^rb 
M  -CnXa 


)  ^-e  , 
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^  ►  X  IT  ^  A  0  f  r  ^ 

^ V  [TSrCv^-  4  wNw  wN ^V,., 

0*/&.  I 


+  HVn  •  W,,^.,*  -W-, 


^N+I^AI 


(Ilk) 


~t'Vi'*'rAM  -  awN  +  ^'^N'rN-rW'Na.|]c^8  +  [f  <W 

*l/j|^|^  |  +•  1*.  r  O  0 


Lksa  & 


le 


hqs .  (113  )ar.d  ( il1*  'are  the  energies  cf  the  lower  edge  half-segcent  before 
the  boundary  conditions  are  applied.  The  integrals  in  (ll^)and  (ll**)can  be 
evaluate!  as  was  dcr.e  previously.  Evaluating  the  integrals  yields  the  following 


energy  expressions: 


-  I -l) 


(V^  * » ^  |  ^  ( a.*  -  a  u.N  V.* 0-~°+  aK^ 

(.1- C&Cy^Ocy*© +-&!/*.  9  AlU\. 

+r»(,40)  M.NWNt  *>~U=  -U.NU.J][(,  -c&& 

+  ~CAc'^ CJ&C  &"  ■ ‘-6  -i**  4«vW-  ^  i=-.]  , 


{115) 


•*£.  -*1  ^'WWN+aw'N+,,^-,+  4l  « 


"H -1*^  ^  r  Ij  |^tl-c«wii^c«<v0+  8  /£*.! 


(116) 


^0  B-h  ft  Ju/A-  +  ^  ^  (.W^  "  V^V-| 

W(M-  /-  B*  Awy\-  B  Aa//\_  '/«*  +  d*\-  t&AX 


-  Q-*\  *t XX,^\_, 


Up  to  this  point,  all  work  on  the  (j=N)  lower  edge  half-segment  has 
beer*  general  and  can  be  used  for  any  edge  condition.  In  order  to  obtain  the 
eqs.  of  motion  for  a  particular  b.  e.,  one  must,  now  specialize  (115)  and 
(llfc)by  use  of  the  appropriate  boundary  condition. 

The  inertial  energy  term  for  the  j=N  point  is  now 

(JV)  \  ,N  *  “  TTf  W  A  LO*  ^  U-^4  1  f-  C|-  oA/A 

^  J  L 

(117) 

Equation  of  Motion  for  the  Point 

The  b.  c.  for  the  guided-pinned  case  are  given  by  (^0)  and  UiU)  will 
be  used  m  eqs.  (115) and  (116) while  (bo)  will  be  saved.  Eqs.  (lib)  is 


t ij  -  ^  (ck-2-)o^t  Ib^-W &J~Wkl  .  • 

(nu-  |  N  M-l 

(118) 

Therefore,  the  bending  energy  for  the  pinned  boundary  edge  half-segment 
is  given  by  equation  (ll6)  i.e. 

|  ?  (o.-Uo(St§)  U.*  f^S-TTlol-iuNfcp) 

u  ~ b) U~N  ^  (p<- 


G*-  ^)lln  >y  o&Lb 

Si}  -  3 

4  >/  G=k-OcjP^-8)t^n_|  -^* 

4-1F  (m  ufN  cefr*  0-H-urN-,lXN^&^  w£  ^  [-aU.  6  & 

'(i - ce^-kG)  9^4  ^  "^p"  (<*- 1)  U&.  5)  U.^  cj&t  8 

-'^zp*  C^“^^r^a)z<JrN_1lJLNosS-5--  L*“^ cA © 

-v "4p-  Cov-  cjf*.  b)  tu*^  -t-1^*  (“A U*J  Ubt  9 

-  T^,y]  [c^b.^  -h  G  -  C^)is)  /k^BJ^P?  (li*  ce*?  & 

UT^)]  [( /-  cb*  &)  ce*  0-  -*VV_  8 

-f  ^>Y\_i“a^^v_  §  -  CL/rv~tQ_/fv_  •  (319) 


Now  noting  eq.  1*0  to  be 


-urN=u_Nti*e; 


(1205 


and  substituting  into  the  energy  expreaaiona  to  obtain  th«  energy  of  the  lower 
edge  half-segment  with  a  guided-pinned  boundary,  the  Bemcraae  energy  tecooea 

®V»  '*?{  KCui -»*. 

^  ■¥  <^C I  f  li)  U*f^  &n  0AtA\ 

(321) 


5^ 


MrC)  G&objl-  fali-fz?)  tart  &  +■  OeK  (  H_  ^ 


v- 


J][ 


fa*  tim  ^  y<*)  ^^.sj 


A.  - 

+  u?  [0-  ^  ^ ^  T^J 


Using  (  iq  )  in  (117)  gives 

)• 


(JO.  *  -  TroWo^  IAjP1 

4 


(.  j  +-  C<*?  £f)  £ajl*v-  0  ~  CW^CJRtibj^# 


P.  (122) 

To  obtain  the  single  eq.  of  motion  for  the  lower  edge  node,  one 


considers  only 


c)v 

olU. 


=  D 


(123) 


N 


There  is  only  one  eq.  of  motion  lor  the  N-node  since  ¥»  +  w^  are  not  in¬ 
dependent  as  shown  by  (  Thus,  one  finds 

(V„V-®-X..N.t.tri  _N  ,(A1  „+(.lrv  ,N_ 


Thus  using  U.21) in  (123/  the  eq.  of  motion  for  the  lower  edge  node  is 


tf-N-l 

(12k) 


given  by 


'W. 


^A)  .  N 

+  - — , 


N 


(125) 


U- 


N 


Substitution  of  the  various  energies  into  (11,) yields  the  equation  cf 
motion  for  the  lower  edge  node  cf  a  guided-pir.ned  boundary  as  follows: 


“■  1  u.  a 

N  w -1  r\i„\  ™  N  N  N  N-a 


(126) 


)vr  -Crr\  ^  W.-  0 
ly  1^-1  N-f  r><  N 


where  . 

^  p  ^  [t  I- 

-*K  z)[u&a  ft  )^M-(j- %*)  5] 4- 

+  9,<*-d  C  Sj'a  w,  /4<|  j 

N  *  /-t)  +  A(i+o)(fe*f  0^<»<+a  SUf  ©)j Jxawv.0 a^>a^k 

f^ipTj-1-  ^J^t3-*va(/*  i>)  To*  $][r,F*>0  & 

A,  ^A>yv  ftjt  <*  jfU-  OSxJ^V^V-f^.-  <U^v_  ft  /^V  y^ 


+  La^H-W  ^  0w\  tftA. 


,  °*  *  '/  \ 
»  A*. - _  r-*  -44M  /OSJ 


.  (  - -  ~  J 

P-ttk'  2L  r  ;-f,  _sA  ,  _  . 

*K*6-o)‘*A  \  ^  ©[C^  ^cje^  ©)  i-ci  <*(?»-  zJ  oe^.0) 

+  os  J  £)  ^-O-cea^V^Sj+O^&r  e^)c^8 

” <-':>^ [c~ P> ^ ■+■ ( i- cpyi  )4. ) ^a^v-  -f  ^  7^>f  ^Pi  £C| 

C^^34>OcfV^  0  -  J^sY^,  &  Aw^,  <^>'t’  JL^v*t5uW-  “g^—  C-AvtA<vv^^OK^ 

,CSZ?  f-  ^rl  ^Lfib'L  ..  r  ^  V  a-n  ^(.N-0  ,  >. 

^  Fjak  ©J^u^vn.  «=*s  _  Sin  jjfc)  "■  Fa+1  PjoS^  ~o^  4*aKi 


+  a  £^i  *  Pj  [a>^  ‘^r^}  xLa  ^  T^v^—  -  WW^ 


CN-l 


A* 


£*  e]c~r-^  , 

*£^7t|  ^*«Citi))Jy^.  >4,  l 

■^  /  ** 

I  ^  £*'  fr[U*-t)c*&S)(5X-  at>  c^*©) 


"<=*(ok  -  CerL  J  ^C  I 

-  c®<*  }U)  (  fv*£  -  P)  't^-“  Zx*  &  jV'- 

-*  ^  C  b)+  &(s*-^u£&L&y  3cxJ  )4* 

“2^ 

1“  C&^i,  /a<}  ZaA.  0^j-  ^  ZS*  0  j^(_  I  ~  006  )/&•')  U&b  ^ 

■  x  & 

&  ^£^vm^.“a  -  J 

-Q^  A*  ^  &-[<£**&, 

S«-  fi  ^S-Avr*.-^  jtas^^c*.  ^CLA-rtcuv\.  ^*5T  "ii-oUi-rt^ferjl 

^\>»  M  '  Trp^'  exT  Ui3  (  !+  6“  S^&OcAO  fiTj  . 


Equations  of  Motion  for  the  j=N-l  Point 


The  eqs.  of  motion  for  the  J=N-1  segment  can  be  found  from 

—  =.  o  . 


(127) 


=  0  . 


(128) 


For  (127)  one  notes  that 


v^.^-eurih  IV. >V+UVN-' 

«  '  (129) 

ind  for  (128^  .  /  rv'N 


Thus  using  (127)  and  (130)  in  (128)and(l27)  respectively, the  eqs.  of 


motion  become 


D(.V- ,) 


o>,A 


N-\ 


v  -  o 

(131) 


(5  tat 


(132) 


Substitution  of  the  various  energies  yields  the  uJJ_1  equation  of  motion 


as  follows : 


CS£,  HHs)  1^-,^  +  (-^N-l  N-|)aN-'  +  ^h-i  n ) 


(133) 
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+t &)[+(-<■- vc&t e)6*  -  s)-°t/w -ai^c*te) 

-  -c*  (_<=*-  z)caL^)+o<  c&^Q 

■^6ava.)^]-^(&  S)[+C«-5v>c^:g)+a^-z^et(9) 
-^][<J^Bu^  /**  ©£(.| 

CQ<i  ^<  ^  /3  ><*t 

-JU\~t&//^ 

i* 


£**/  ©^s  flS/r\,  4< 


£Lc* 


^ y^ggi? 

"ZTvl^A-  ^  %K+Y  F(j* 

s<  ^ 

-pt)c^.e)  -a«C=*-ai)c©te)][^5  msk.Z< -Q 

-  ce*>  K*)  bJ+z) 


*3 


C^<) 


^N-ll  , 


^*vyv-w  f.  2z)  e*  006  / 


e»<  /ck 


t*  ^T 


t 

+  ()-«»)£OxOw  &] 

-  Wfccu^  slr]+ff6c00<?)j-a 

A. 


+  JUCW  ^r-(Ut^^i]+^[o-ca6>^)x 

0  ~  h  t 

and  "M 

p-^W  _  .  a  A  •  5lCN-|)  |  y 

'  N'l  N'l  ""  ^  77^  rV.C>w  U')  ^iAyTw  ex  T^L/n^  j 


Equations  of  Motion  for  the  j=N-2  Point 

At  first,  it  appears  that  the  2nd  point  off  the  boundary  of  the  shell 

would  be  a  general  point.  However,  as  in  the  case  of  the  dome  point,  this 

is  not  true  of  the  j=N-2  node.  The  eqs.  of  motion  of  the  lower  edge  node 

u  w 

except  for  the  coefficient  Sjj  ^  g ,  since  all  coefficients  except  this  one 
check  out  as  symmetric,  are  given  by  the  general  eqs.  of  motion,  i.e.  eq:_ . 
(108).  Thus  the  only  eq.  to  be  formulated  in  this  section  is 


-&LY. 


0 


-0 


vj 


(135) 


For  (135)  one  can  note  that 


to  I ^M-a+ 


I 


+<JtV 

Thus  using  (136) in  (l37)gives  the  equation  of  motion  as 

W, ..  5cvav  uo  £>cw»>,  „ 


(136) 

3(Jl) 


rO 


cHv, 


d 


N-3.  "N‘J>  v-,,'t‘rN-»  — n-i 

Evaluation  of  the  various  coefficients  defined  yields, 


-OOT 


d 


(137) 


/  -  -U/"  ^  \  /  -OU"  tA-  \  UA-\  -W  \ 

^ ^  N-a’  N-r  N +  N-«' 

x  n/J~N_4+  ^  N-a  iv- 3)  N-a^^v-a1"  ^  N-a  n-  i 


^  ^  N  *  3v  N>»  ^”{\j  -  9s,  ” 


D 


(138) 


where  l. 

-w  uc 

^  Al  ^  M _ ^  ^  / 


Ai-JK  N-"3l 


X/u^ 

-iaT  u.  KfA^TT 


c^-r  ,  .  atN-al  .  .  v 

“0  0&4)(sX  tl  +  u)c&> 


SiLN-z) 


^K-JN.,  3  A 0  |_a»Cl  +  T>),A^.  ex  C&l'/U J 
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-W  U.  1  ^  .  \y 

N  “  |"V  '  72rt  B)ju^  &<  A*s^  jC< 


+  ( fcn  .&)  o^yv-i)  y^< 


cw  -w  r^*1  ,  .  x 

■^N-Sv  N*4  "  /^(/-0)aT I  <y  Ju*^,  <=*  Aa/a-A 
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PROBLEM  FORMULATION 


Summary  of  Equations  of  Motion 

Since  all  the  required  equations  of  motion  have  now  been  derived  it  is 
possible  to  gather  all  the  eqs.  of  motion  for  the  guided-pinned  spherical 
cap  for  easier  reference.  They  are  as  follows: 


je0:  eq.  of  motion  w.r.t.  w  (dome  point): 


csr  ~s)  1-  cs^~)  --  ttsr  r)^-Lvcrhr°> 


J=1 :  eq.  of  motion  w.r.t.  u  : 


(139) 


eqs.  of  motion  w.r.t.  w 


(lbO) 


(sfr)  ,4(  $  r  & r  ?)  fcrrv,  ^ ; 

-  aw 

J=2 :  eq.  of  motion  w.r.t.  u  (with  k=2):  * 

« 1 7)  K.,  Ks  3-  3:  )  «_AK  ss  S')  -3+  CS  £  'f)w,  +  Cs  r X ^ 

)  (H2) 

eqs.  of  motion  w.r.t.  w^: 

tssr  t )  «.  sr  t)  + »  £■  sr  H+(  ^  ns 11131 

t  nr  1  Li^1 x )  ■ t  *n*r‘i- -wz  *)*$=  ot 

J=k:  (k  =  3jii . «-**,  N-3): 

eqs.  of  motion  w.r.t.  u 

St,.,') (.^KV.W^-r  fc-)ur.lk-i 

eqs.  of  motion  w.r.t.  w,  : 

(^X  h)  ^>+CSX  £r)u 

>T -f-Cs.K  CsX 

JgN-2:  (11*5) 


(lM) 

,-ur-ur  \ 


eqs.  of  motion  w.r.t.  u^  ^  wilh  k*N-2: 
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eqs.  of  motion  w.r.t.  w 


(C/^  ul  \  ,  /-w  ^  \  u.\ 

Ai-a  N-3^V aj-* 
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eqs.  of  motion  w.r.t.  u 
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Matrix  Representation 

It  is  evident  that  most  of  ...e  _t:ff..e„s  coefficients  have  a  symmetric 
property.  Hence  the  eqs.  of  motion  can  be  arranged  in  such  an  order  that  they 
will  have  a  symmetric  stiffness  matrix.  From  McDonald's  paper  [yj  one 
observes  that  there  appears  to  be  two  methods  of  arrangement  for  the  dis¬ 


placement  vector.  They  are  for  our  case 
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The  first  type  looks  the  most  promising  since  to  obtain  the  second  type 


vector^one  would  have  to  put  the  3  eqs.  of  J-0  and  J=1  in  some  symmetrical 
form  and  this  seems  to  be  impossible.  Therefore,  consider  all  the  u  eqs. 
and  then  all  the  v  equations.  Hence,  one  obtains  the  following  stiffness 
matrix: 
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Therefore,  the  system  of  2(N+l)-2  eqs.  of  motion  can  be  written  as 

[Cs]  -O]]  [%]  ■ 


(151) 


How ,  one  can  group  into  dimensionless  quantities  and  rearrange  terms 
in  eqs .  ( i 5 ]. )by  multiplying  them  by 

(.  I  ~  ill -  . 

h~7TK  J 

so  that  the  group  coefficients  of  the  energy  expression  are  as  follows: 

''.‘‘zH--  H-  ■  / 

:  faiM 3a?  I  '  '  1  I??'  . I  . 


Hence  eqs.  (1)  car.  be  written  as 


)c- ]-*'[>] I  cti*0 


where  [cj  is  now  the  stiffness  matrix, 

1  i  A  A 


X1  -- (f|  7(\^)(X-X), 


(153) 


(151*) 


where 


*.  E/a  ( i+£) 


Therefore  the  new  mass  matrix  is  of  the  form: 


Now  define 
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au-o*) 
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[A]  • 


(156) 


f 

Now  the  stiffness  coefficients  of  the  stiffness  matrix  from  the  eqs.  of 
Dtion  will  be  redefined  to  coincide  with  eqs.  (153)  Thus  the  coefficients 


are  listed  as  follows: 
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The  coefficients  of  the  u  eqs.  are  given  by  the  general  expressions  for  k=2,3..., 
N-3,  N-2,  so  that 
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